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obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  fommilated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
<  .<her  person  or  corporation,  or  conveying  any  rif^ts 
ui  pcinissinn  to  manufacture,  use  jr  >5011  any 
patented  invention  that  may  in  any  way  b?  i-ele+ci 
thereto. 
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THEORY  AND  DESIGN  CURVES  FOR 
A  YO-YO  DE-SPIN  MECHANISM 
FOR  SATELLITES 

by 

J.  V.  Fedor 

Goddard  Space  Plight  Center 


SUMMARY 

Equations  are  developed  fur  the  yo-yo  de-spin  mechanism. 
Under  conditions  usually  met  in  a  de-spin  applicofioa,  ihc  equa¬ 
tions  can  be  greatly  simplified.  A  computatioi;  sheet  is  pi  atad 
based  on  Ihe  simplified  equations.  The  design  of  a  de-spin  mecha¬ 
nism  is  thus  reduced  to  reference  to  a  chart  and  a  few  simple 
calculations. 
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THEORY  AND  DESIGN  CURVES  FOR 
A  YO-YO  DE-SPIN  MECHANISM 
FOR  SATELLITES 


by 

J.  V.  Fedor 

Goddard  Space  Plight  Center 


INTRODUCTION 

The  yo-yo  de-spin  mechanism  is  essentially  two  pieces  of  wire  with  weights  on  the 
ends  (Figure  1).  These  wires  are  symmetrically  wrapped  around  the  equator  of  the 
satellite  and  the  weights  are  secured  by  a  release  mechanism  (Figure  2).  At  a  pre¬ 
selected  time  after  satellite  spin-up  and  separation  from  the  launching  vehicle,  the  weights 
are  released,  thue  discarding  enough  momentum  to  reduce  the  spin  of  the  satellite  to  the 
desired  value.  It  is  Uie  object  cf  this  note  to  pccscnt  tha  jo  .-c  di— sp’n  theory  in  eqnat’cn 
and  curve  form  so  that  the  designer  can  readily  apply  the  theory  to  e  dtsi^n  of  {no  de-spin 
msehanisr  , 


Fc-  .1  values  of  the  satellite  moment  of  inertia  and  radius,  c-nd  for  a  given  amount 
of  spin  reduction,  the  designer  wants  to  determine  tiie  weights  ihai  i<uot  t'*"  use  .1,  'Jio 
length  of  wire,  and  the  maximum  tension  in  the  wire.  In  what  follows,  equations  will  be 
developed  for  these  quantities.  It  will  also  be  seen  that  under  certain  conditions  usually 
met  in  de-spin  applications,  a  great  simplification  of  the  equations  is  possible.  With 
this  simplification,  all  cf  the  design  curves  can  be  put  onto  one  dimensionless  graph. 


ANALYSIS 

After  being  placed  in  orbit,  the  satellite  travels  at  essentially  constant  velocity  and 
the  centrifugal  action  is  cancelled  out  by  the  gravity  force.  Since  any  equation  of  motion 
written  for  the  satellite  would  be  independent  of  this  constant  velocity,  the  satellite  can  be 
considered  stationary  and  spinning  about  a  fixed  axis  for  the  purposes  of  the  present 
analysis. 
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There  are  two  phases  to  the  spin  reduction  process  considered.  In  Phase  1  the  wire 
is  changing  in  length  and  is  tangent  to  the  satellite.  In  Phase  2  the  length  of  the  wire  is 
constant  but  its  position  is  changing  from  tangent  to  perpendicular  to  the  satellite.  With 
proper  design,  the  wire  can  then  be  released  when  it  is  perpendicular  to  the  satellite.  It 
should  be  noted  that  releasing  the  wires  after  Phase  2  is  more  efficient,  weight-wise, 
than  releasing  them  after  Phase  1. 


Phase  1  Analysis 


A  sketch  of  Phase  1  with  its  coor  Jmate  system  is  shown  in  Figure  3.  Only  one  wire 
is  shown,  as  tt.e  system  is  considered  symmetrical.  Also,  the  system  is  considered 
torque-free;  small  moments  due  to  the  earth's  magnetic  or  gravitational  fields  are 
neglected. 


7.'he  total  kinetic  energy  of  the  system,  and  also  the  Lagrangian,  is 


T  =  ^  +  I  in(x2  +  y^)  , 


(1) 


where  is  the  'otai  mass  of  both  weights  and  t  is  the  m''.nent  of  inertia  of  the  satellite 
about  the  spin  axis.  Initially  the  wires  are  considered  weightless;  Jacer.  ?  -.rr-fiu/u  wil. 
be  d'.veioped  to  fcu;e  the  weight  of  the  wire  into  account  without  disturbing  the  form  of  the 
Gouat’f'ns.  U.-j’.ng  the  transformation  equations 

(2) 
(3) 


X  =  a  cos  6  i  si:;  ^  . 
y  -  s  sin  9  ~  t  cos  9  , 


the  Lagrangian  can  be  put  into  the  form 

T  =  -I  10^  +  .  (4) 

The  length  of  wire  unwound  at  any  time 
during  this  phase  is 

t  =  a(9  -  <t>)  .  (5) 

The  equations  of  motion  in  Lagrangian 
notation  are 


d 

dt\?i 

[\  _  31 
hj  30 

0 

II 

II 

(6) 

dt\^ 

i\  3T 
?/  ■  30 

=  %  =  0. 

(7) 

y 


r  oure  3  -  Phase  1  coordinate  system 
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where  <f>  and  e  are  considered  to  be  generalized  coordinates.  Equations  5,  6,  and  7 
can  be  integrated  once  to  yield  an  equation  which  e:q)resses  the  conservation  of  momen¬ 
tum.  Equation  4,  v/hich  is  equal  to  a  constant  and  represents  the  conservation  of  energy, 
can  be  combined  with  the  conservation  of  momentum  equation  to  yield  ^  and  0  as  a 
function  of  the  length  of  the  wire  during  Phase  1: 


.  0„(i  - 

*  =  “TTTW 

;uid 

24 


(8a) 


(8b) 


where  4  is  the  initial  spin  rate  and  =  i  /m  +  .  The  above-mentioned  equations 

can  also  be  integrated  to  give  the  length  of  wire  as  a  function  of  time: 

«(t)  =  a4t  .  (9) 

The  total  tension  in  the  wires  is  simply  the  product  of  the  mass  of  the  weights  and  the 
uccoleration: 


”  -m(a<?  +  £^2)  ,  (iOa) 

•,Tiie."C  'S  the  acceleration  of  the  mass  m  along  the  wire.  Substituting  Equations  8a  and 
8b  into  Equation  !0a  yields 


^  4m4^£  (1  -  a^/\2) 

"  (1  +  e2A2)='  ■ 

It  can  easily  be  shown  that  the  maximum  tension  during  Phase  1  is 

Fi  =  |V3mX42(l  -  aV.'.2)  , 


(lOb) 


(11) 


or  in  dimensionless  form 


mcoQ  \ 


1.3(1  -  aW) 


where  =  4  • 


(12) 
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Phase  2  Analysis 

A  sketch  of  the  Phase  2  arrangement  is  shown  in  Figure  4,  The  total  kinetic  energy 
in  this  phase  is 

T  =  -i  1.^2  ^  1  in[a202  ^  jat  cos  {6  -  y)  4>y  +  t^y^l- 

where  ^  ^  (from  Equation  5)  since  t  is  a  constant.  This  phase  is  much  more  com¬ 

plicated  analytically  than  the  previous  one.  Fortunately,  we  are  mainly  interested  in  the 
spin  at  the  end  of  Phase  2  (at  the  release  of  the  wire,  when  ^  -  y  =  0)  and  ti.e  •.  .;sion  in 
the  wii'.%  They  beth  can  be  obtained  by  applying  the  conservation  of  momentum  and  of 
energy,  and  an  equation  similar  tr  iSquation  10. 

The  conservation  of  energy  (at  0  -  y  =  o)  gives 

^  j  m(a^2  +  ty)2  =  const.  =  i(I+,na2)<^2  (14) 


and  the  conservation  of  momentum  gives 

i  !5(f  +  (3!/>2  T  Cy)  -  const.  =  (I  +  ma2)(^jj. 

Combining  Equations  14  and  15  we  obtain 


(1  -  r)I 


+  1 


y 


'(1  , , 


(15) 


(16) 
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In  terms  of  g  ,  r  and  e  /a  ,  the  force  is 


^  9«/ar2  Tg  H  1  -  r(t/a  + 

Fj  -  +  l)e/a  J  J  ‘ 


The  force  given  by  Equation  18  or  19  is  a  maximum  for  Phase  2  because  the  length  of  the 
wire  is  constant  and  the  angular  velocity  of  the  wire  (y)  is  a  maximum  just  before  release 
(the  a  r2/f  term  contributes  very  little  to  the  force). 


Simplification  of  Equoiions 

The  equations  developed  thus  far  are  exact  in  that  no  analytical  approximations  have 
been  made  in  their  derivations.  A  great  simplification  can  be  made  by  noting  the  follow¬ 
ing;  G  in  Equation  17  is  large  (in  the  order  of  200  or  more)  for  most  practical  spin  weights 
and  spin  reductions.  Thus  if  G  is  replaced  by  G  + 1  in  Uie  denominator  of  Equation  17,  we 
have 

-  (G .  1) 

(t/a  ,  i)-  G  +  1  ,  (20) 

r  s,  .  I  ——I  -  s  ■■  -  1  . 

G  -  1  (1'/.^  +  !/ 

Re?;  raiiginr  Equ^tior.  20  into  a  dimensionless  form  gives 

.  J_-  r  1  t  r  -  (g/a  t  1)'^  ,21) 

ra(a  +  E  1  -  T 

Now,  if  E/a  is  greater  than  2n  (t  is  one  circumference  length)  the  (E/a  +  1)“^  term  in 
Equation  21  is  negligible.  Hence,  we  have  the  very  simple  expression 


in(a  +  E  1  -  r  ^ 

Equation  22  has  been  plotted  in  Figure  5. 

Calculations  indicate  that  the  maximum  tension  in  the  wire  occurs  in  Phase  1  in  most 
cases;  hence  the  simple  expression  of  Equation  12  can  be  used  to  calculate  this  maximum 
force.  Since  a2/\2is  small  in  most  de-spin  designs,  the  maximum  force  is  simply 

=  1.3mc.2\  ■  (23) 

Consistent  with  Equation  22,  Equation  19  can  be  simplified  to 

far2  [(r  +l)(l/a  +  1) 

P_  =  "''■^0^  —  +  ! - - 


Distributed  Mass  of  Wire 


We  shall  now  develop  an  approxi¬ 
mate  method  to  include  the  distributed 
weight  of  the  wires  without  disturbing 
the  form  of  the  derived  equations.  From 
Figure  6  tiis  kinetic  energy  of  i***  par¬ 
ticle  is  i  (xf  +  y?)  ,  where  p  is 
the  mass  density  per  unit  length  of  t’  e 
wire.  Sumr  ’ng  over  all  the  particles, 
we  ha/e 

i  X]  pAi-  (xf  +  yf)  .  (25) 

By  using  the  transformation  equations 


Xj  =  a  cos  5  f  cos  y  ,  (26) 

y.  =  a  sin  +  f  sin  y  .  (27) 


Figure  6  -  Coordinate  system  for  distributed 
weight  of  wire 


and  replacing  the  summation  by  integration,  the  total  kinetic  energy  of  the  wire  is  found 
to  be 


®wire  =  I  r 

'*  *0  ' 

=  i  p^fa^^^  +  af^  cos  -  y)  4^y  • 


Since  t  is  large  just  before  release,  —7'-  is  rhe  large.st  term  in  Equation  10.  in  any  ap¬ 
proximate  formula  we  should  include  all  of  this  term.  Separating  Equation  28  into  two 
parts  results  in  the  following: 


®wire  =  '5  cos(e  -  7)07  +f^7*J 

+  ^  p|-  +  X  cos  (6  -  7)  07]  • 


Just  before  wire  release,  0  is  also  small,  so  the  seccr.-’  part  of  Equation  29  is  small 
compared  to  the  first  part.  Hence, 

*®wirc  ~  I  T  ^  cos(e-v)07 

where  p£  /3  is  one-third  the  mass  of  the  wire.  Note  that  Equation  30  is  similar  to  the 
second  term  of  Equation  13.  Note  also  that  when  B-y  =  v/2  and  y  =  e  (that  is,  Phase  1), 
Equation  30  is  similar  to  Equation  4.  With  this  approximation  in  Phase  1,  the  kinetic 
energy  due  to  the  change  of  wire  length  is  ai  ?o  neglected.  The  change  of  length  of  the 
wire  does  not  contribute  much  to  the  total  kinetic  energy  of  the  wire.  Thus  we  can 
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include  the  effect  of  the  weight  (mass)  of  the  wire  by  including  one-third  the  weight  of 
the  wire  in  the  weight  of  the  mass  term  of  the  derived  equations.  We  therefore  define 
m  in  the  derived  equations  in  the  following  way: 

pt 

m  =  nio  + 

where  niQ  is  the  mass  of  the  spin  weights  and  p  is  twice  the  normal  mass  density  of  the 
wire  (to  take  into  account  both  wires).  This  equation  should  be  valid  for  large  spin 
reductions  (r  small). 


APPLICAiiON  OF  DESIGN  CURVES  AND  EQUATIONS 

We  shall  now  apply  the  simplified  equations  to  a  practical  case.  In  Reference  1,  the 
following  information  is  given  for  the  S-30  Ionosphere  Probe  Satellite: 

I  =  2.5  slug  ft^ 
a  =  1.25  ft 
I  15w  rad/sec 
"f  =  10?r/3  rad/sec 

A  calcuhuion  for  r  gives 

I  r  10/45  =  0.222 


From  Figure  v,  with  r  =  0.222,  we  read  the  value  i  57  for  i/m(E  +  a)^.  If  we  let 
t  s  '’00  ineP.'?.'.  a  .simple  calculation  for  the  weight  of  the  spin  weightr  md  wire  gives 


w  .:  ing 


Ig 

1.57(e  +  a) 


1 


2.5(32.2) 


0.151  lb. 


A  check  of  the  assumption  that  G  is  large  give's  a  value  of  266  which  is  more  than  ade¬ 
quate.  Reference  1  reports  a  value  of  0.150  lb.  calculated  from  the  theoretically  correct 
equations. 


Taking  into  account  the  weight  of  the  wire,  we  have  for  the  spin  .v-^ights  alone 
niog  =  '"R  - 

=  0.151  -  i  (0.005)  =  0.1211b. 

where  a  double  weight  density  of  0.005  lb  per  foot  has  been  assumed. 


The  calculated  maximum  tension  in  one  wire  is 


=  1.3 


1.3(0.151)  (47)2  (23.12) 
?(.32.2) 


156.6  lbs 
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The  value  read  from  a  graph  in  Reference  1  is  155  lbs.  A  calculation  of  the  maximum 
force  in  Pliase  2  (Equation  24)  gives  150  lbs. 

To  aid  in  systematizing  de-spin  calculations,  a  computation  sheet  has  been  developed 
(Appendix  A).  It  has  been  found  convenient  to  use  in  de-spin  calculations.  It  is  believed, 
also,  that  this  sheet  helps  to  reduce  human  errors. 


CONCLUSION 

Th  sirjplified  equations  derived  here  greatiy  facilitate  the  application  of  :;c-spin 
mechanism  theory  to  a  particular  desig.-u  With  the  calculation  sheet,  the  procedure 
has  been  reduced  to  reading  a  graph  and  making  a  few  routine  calculations. 
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Appendix  A 

Yo-Yo  De-Spin  Cciculation  Sheet  (Radial  Release) 


DEFmiTION  OF  SYMBOLS  AND  UNITS: 

I  -  moment  of  inertia  about  spin  axis  (slug  ft^) 
a  -  radius  of  satellite  (ft) 
e  -  length  of  one  yo-yo  wire  (ft) 

m  -  total  mass  of  both  spin  weights  +  1/3  mass  of  both  wires  (slugs) 

F^.  -  maximum  tension  in  wive  (lb) 
cj'  -  initial  f:p{ri  rate  (vsn\/sec) 

-  final  spin  rate  (rad/sec) 
r  -  final  spin  rate  divided  initial  spin  rate 
g  -  acceleration  of  gravity  (ft/sec^) 

TO  CALCULATE  THE  TOTAL  MASS  (WEIGHT)  OF  SPIN  WEIGHTS  AND  WIRE  (m): 
Record 

I  = _ slug-ft^  =  _ rad/sec 

a  = _ ft  ^  =  _ rad/sec 

Calculate 


With  this  TOlue  of  r  ,  read  the  value  of  I/fe(£  +  a)*  from  the  design  curve;  <.  ill  this  value 
B.  Then  calculate  the  followina: 


B.  Tiien  calculate  the  following: 

ig 


'.V  =  mg  = 


=  J _ >32.2 

B(e  +  a)2  (  H  ^ 


lbs. 


TO  CALCULATE  MAXIMUM  TENSION  IN  ONE  WISE.  Calculate  K  by 

=-L+  a2  =  _  =  _ . 

m 


or 


Also 


ft. 


„2  = 


/sec^: 


1.3  =  1.3  ( - ^ 


)=  ( 


)  = 


lbs. 


CHECK  OF  UNDERLYING  ASSUMPTION  OF  THE  EQUATIONS:  Calculate  G  as  follows: 


G  = 


If  G  >  100  and  f /a  >  277  ,  the  answers  are  accurate  to  about  1-1/2  percent  of  the  theo- 
reticaUy  correct  value. 
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